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Abstract. This paper proves that if points Z\,Zi,.., are chosen independently 
and identically using some measure /z from the unit circle in the complex plane, 
with p n (z) = (z — Z\)(z ~ Z-i)...(z — Z n ), then the empirical distribution of the 
critical points of p n converges weakly to fi. 



1. Introduction 

Across many fields of mathematics, one of the fundamental questions about a func- 
tion is the location of its zeros. Perhaps the most famous example is the Riemann 
Hypothesis. Beyond that, entire fields such as algebraic geometry and the emergent 
study of stable functions have locations of zeros as their focus. 

The relation between the zeros of a function and the zeros of its derivative (the 
critical points) is interesting and not always obvious. In the case where all zeros are 
real, Rolle's theorem tells us that the zeros of the derivative interlace the zeros of 
the function itself. In the case of complex polynomials the analogous result is the 
Gauss-Lucas theorem which states that the zeros of the derivative of / must lie in the 
convex hull of the zeros of / and gives a representation of the zeros of /' as convex 
combinations of the zeros of /. A corollary of this is that differentiating preserves 
stability. Differentiation is also known never to increase the number of non-real zeros 
of a polynomial. 

Two famous conjectures in this area are the conjectures of Sendov and Smale. 
The former, made by Blagovest Sendov during the 1950's, states that if the roots 
Z\, Z2, z n of a polynomial all lie inside the closed unit disc, then for each root of 
the polynomial, the closed unit disc centered at the root must contain at least one 
critical point. The latter, made by Steve Smale, states that if / is a polynomial of 
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degree n with at least one root and /' 



(0) ^ 0, then, 
: /'(£) = o) < K, 




where K — 1 or Sendov's conjecture has been proven for the case when 

zi, z 2 , z n all lie on the unit circle, whereas Smale's conjecture has been proven 
for when / has all its roots, save 0, on the unit circle. The most general forms 
of these conjectures are still unsolved. More information on these conjectures and 
proofs of some of the special cases can be found in [RS02J. 

A probabilistic study on the roots of derivatives of polynomials was done by Pe- 
mantle and Rivin in |PR12j . Let / be a polynomial with n roots that are chosen 
independently and uniformly from a measure [i on the complex plane. They con- 
jectured that the empirical distribution of the roots of /' converges weakly to \i as 
n — > oo. They prove this in the special case when \x has finite 1-energy, namely when 
\x satisfies 



This condition cannot hold, however, when \i is supported on any set of dimension 
1 or less. The aim of the present paper is to extend their result to the case of any 
measure supported on the unit circle. 



Say, Zy^Z-i, ... is a sequence of points chosen i.i.d. with respect to some distri- 
bution fi on the unit circle. Write, Z^, = exp(27ri6k), so that is a collection of 
IID random variables whose common law is supported on [0, 1], which we denote by v. 



and y{ n> , yf>, y£i x be the roots of p' n (z). 

For k > 1, let = E(Z k ), where Z ~ fi. We know, by the Strong Law of Large 




2. Notations and Background 



Let 



p n (z) = {z- Zx)(z - Z 2 )...{z - Z n ), 



Numbers, 




Z k 1 +Z^ + ...Z k 



a.s. 



n 
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Denote by Z(f) the empirical distribution of the roots of a random polynomial /. 
That is, if / has roots X 1 , X 2 , X m , then Z(f) = ~ Y!j=i &Xj- 

We shall write D for the open unit disc, and C for the unit circle. 

In their paper, |PR12j . the authors conjectured that, for any distribution \i on the 
closed unit disc, Z(p' n ) converges weakly to p. That paper also proves the following 
proposition. 

Proposition 2.1. Let p be the uniform measure on C. Then Z(p' n ) converges to C 
in probability, that is, P(Z(S) > e) — > 0) for any e > and any closed set SCD, 
disjoint from C . □ 

In this note, we shall generalize this to prove that 

Lemma 2.2. For any distribution p on C, Z(p' n ) converges to C in probability. In 
fact, if p is not uniform on C, the convergence is almost everywhere. 

The above leads us to prove our main result, which is a special case of the afore- 
mentioned conjecture in |PR12j : 

Theorem 2.3. For any distribution p on C, Z(p' n ) converges weakly to p on C. 

The proof, as shall be seen in forthcoming sections, can be divided in to two parts, 
the latter following a pattern similar to the proof of Weyl's equidistribution criterion 
(see, for example |Ch68j ). The former requires the following theorem proved in 
[CN06| regarding a companion matrix of the critical points. 

Proposition. If zi, z n G C, and yi, y^, y n -i ar ^ the critical points of the 
polynomial p n (z) — (z — Z\)(z — Zi)...(z — z n ), then, the matrix 

(i) D (i-£\ + *±j 

\ n J n 

has yi,y2, Vn-i as its eigenvalues, where D = diag(z±, z 2 , z n _i), I is the identity 
matrix of order n — 1 and J is the {n — 1) x {n — 1) matrix of all entries 1. □ 



3. Proofs of Lemma 12.21 and Theorem 12.31 

We first begin by proving a small lemma. 

Lemma 3.1. Let p, be a distribution on the unit circle C with Ck = E(Z fc ) ; where 
Z ~ p. Then Ck = for all k > 1 if and only if p is uniform on C. 
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Proof. Clearly if /i is uniform on C then Ck = for all A; > 1. Now say fi is not 
uniform on the circle but we still have Ck = for all k > 1. Then the law v is not 
uniform on [0, 1]. Now, if Z-y, Z2, ... are points on C, chosen i.i.d. using /x, and if we 
write Zj = exp(27ii9j), j = 1, 2, then 9%, #2, ... are points in [0, 1] that are i.i.d. v. 

As mentioned earlier, by the Strong Law of Large Numbers, for all k > 1, 

z k l+ z k 2 +...z k n a . 



0. 



n 

and so by Weyl's criterion, for any < a < b < 1, 

YTj=l 1 {e J e[a,b]} a. 



b — a. 



n 



But l{0 i6 [ Oj &]}, j = 1, 2, ... are i.i.d. random variables taking values or 1 with expec- 
tation ^([a, &]). So by Strong Law of Large Numbers, 



YJj=l 1 {e 3 e[a,b]} a.s. 



V 



n 



(Ml). 



Since, v is not uniform on [0, 1] we have arrived at a contradiction. Therefore, there 
must exist at least one non-zero c^. □ 

We proceed to use this fact for the proof of Lemma 12. 2[ 



Proof of Lemma Iff.M Assume fi is not the uniform distribution on the circle (as the 
uniform case has been taken care of in [PR12J). Then, as mentioned above, there is 
at least one non-zero Cfc. Thus the power series function f(z) = Y^T=o^k+iz k exists 
at every point z £ D, is analytic there (since \ck\ < l,Vfc), and so has only finitely 
many zeros inside any r-ball, where r < 1. 

Define 

np n (z) - Zj 

V n has n — 1 zeros, which are exactly the zeros of p' n (z), and n poles, which are 
exactly the zeros of p n (z). Thus V n (z) is analytic inside D. We shall show that as 
n — > 00, V n converges inside the disc to — /, uniformly over compact sets. To see 
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this, note that for zGD, 



j=l I 3 j=i 3 \ 

1 n oo oo 



n 

3=1 fc=0 fc=0 



where, we write a^ +1 for the fcth power sum average z i +z 2+---+ z n _ ^ g s t a ted earlier, 
a£ Cfc for all fc > 1. 

Let < r < 1. Given any 6 > 0, > 1 such that 



oo 



< 



r 4 

Corresponding to the chosen K, there exists an N > 1 such that, 

i k , 5(1 -r) 

K-Cfc| < 2 - , 

Vn> N and VA; = 1, 2, if - 1. Therefore, Vn > N and all 2 G £ r (0), 

A'-l 00 



|K(z) + /(^)|<^|<-c fc |r fe + £ 



fe=0 k=K+l 



<%-^--(l+r + r 2 + ... + r K ~ l ) + 2. 6 -<5, 
- 2 v 7 4 

which proves uniform convergence of V n to — / over compact sets. 

Hurwitz's theorem (see, for example, [Co78]) says that, if g n is a sequence of ana- 
lytic functions that converges uniformly on compact subsets of an open set G in the 
complex plane to an analytic function g, then, for any closed disc D contained in G, 
if g has no zeros on the boundary of D, there shall exist a natural number N for 
which g n and g have the same number of zeros in the interior of D, for all n > N. 



Using Hurwitz's theorem, given any < r < 1, there exists an M > 1 for which 
V n and / have the same number of zeros inside B r (0) for all n > M. That is, p' n and 
/ shall have the same number of zeros inside B r (0) for all n > M. But, as discussed 
above, / has only finitely many zeros inside B r (0). Thus Z(p' n ) converges to the unit 
circle almost surely. □ 
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Our main result, Theorem 12. 3 \ will be a consequence of the following proposition 
together with Lemma 13.31 below. 

Proposition 3.2. Given any sequence of points z\,Z2,--- with \z n \ < M for all n, 
and z i +z z+--- z ™ — t Ck as n oc, Vk > 1, the critical points y± ,y% > •••;2/i_i °f 
p n (z) — (z — Zi)(z — z-i)...{z — z n ) also satisfy 



(l/J B) )* + (l/J ,) )* + ... + (y& fc 



n 



Ck as n — > 00, 



VA; > 1. 



Proof. Note that, it is easy to see that this theorem holds true for k = I, because 
the average of the critical points is exactly equal to the average of the roots (by 
comparing the coefficients of z n ~ l in p n (z) with z n ~ 2 of p' n (z)). To prove the result 
for general k, we use a result of [ CN06] . mentioned as a proposition in Section 2, to see 
that for k > 2, (y^ ] ) k , (y^ ] ) k , {y ( nli) k are the eigenvalues of [D (I - i j) + f J] k , 
and so, 



{y?) k + (y?) k + - + iv { :\) k = Tr 



D I - -J + —J 



n 



n 



Note that the expansion of [D (J — - J) + — J] k is the sum of all terms such as 



(2) D h 



DJ\ 



n 



1 



— J 

n 



D 



DJ 



n 



— J 

n 



...D 



n 



n 



where the exponents li, 1%, are non-zero integers, with hj-2 + hj-i + hj — 1 
for all j = 1, 2, .., k. Clearly the number of such terms is 3 k , which does not depend 
on n, and so, if we find that the trace of the matrix in the expression ([2]) converges 
as n — >■ 00 to aj 1 ,i 2) ...j 3fc , then the trace of [D (J — ^ J) + ^J] k converges to a h,h,-,hk- 



Henceforth, we fix l\,l2, ■■■hk- Now, note that J r ' 
and 

... z{\ 
... z\ 

■■■ <J 



(D p J)(D q J) 




for any p,q > 1. 



n _ \yn Ijm 1 f or an y m > 1 ; 



/ 4 



4 ■■ 








z\ •• 


■ 4 




' 4-1/ 
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The above tells us that there exists p, q, s , s\, s 2 , s^-i > such that, term ([2]) 
is of the form 

(3) 

where the numbers p, q, s , si, Sk-i are determined solely by the ij's (and so, are 
independent of n). For instance, 

P = h + k + ■■■ + hk-i, and 
q = h + k + ■■■ + hk- 

Also, M can only be one of the following terms: D k or ^— ^ or ° mi J D m2 for some 
m, mi,r/i2 > 0, which are fixed, < fc, and dependent only on the Zj's (this fact, and 
also the scalar coefficient in ([3]), can be checked by simply using induction on k). 
Furthermore, the scalar coefficient in is always 0(1). 

Observe that, if M = D k , then the scalar coefficient of M is equal to 1 (since there 
is only one D k term in the expansion of [D (i — - j) + ^J] k ) and TliMl Cfc# 

On the other hand, if M — then 

rp (M , {y?) m + {yf ) ) m + - + {y i nli) m , , 

Tr(M) = = o{n), 

and if M= ^U-D m \ 

n 

Tr(M) =Tr( D mi+m2 - 
V n 



^( n )jmi+m2 _|_ ^( n )^mi+m2 _|_ _|_ fy( n ) ^mi+m2 



o\n) 



n 



Thus, 



Tr [D (I - i J) +fj] k 

— - — > c k as n — > oo. 



n 



□ 



We require now a simple lemma about the convergence of weighted averages with 
some suitable conditions on the weights. 

Lemma 3.3. Let {a?nfc}njfe=i be a double sequence of real numbers or real-valued 
random variables, for which, x n k G (0, 1) for all n and for all k. Also, let {b n k} be 
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a double sequence of complex-valued random variables, with \b n k\ < M for all n and 
for all k. Then, if 

Xln + x 2n + ■■■ + 



n 

in probability, almost surely, or pointwise, 

x lnbln + x 2n"2n + ••• + X nn b 



-+0, 



n 

respectively in probability, almost surely, or pointwise. 



" n -> 0, 



Proof. We have, 



t" • • • ~\~ x nn b ri 



n 



\- ... -\- x nn o n 



< M 



n 

X\ n "T %2n i ••• "T *Enn 



n 

which gives us the result. □ 
We now have all the tools required to prove our main result, namely Theorem 12.31 

Proof of Theorem \2.3l Say we write, 

Vj = ^ ex P( 27T ^j n) )'i = !' 2 ' 
The proof will consist of three major segments. Our first task is to prove that 

n-l 



n — 1 z — ' J 



In fact, unless /i is uniform on the circle, we will show that 

n-l 



n — 1 z — ' 



Next, we shall use the above information and Lemma [3.31 to show that 
exp(2kni(f) ( ( 1 ^) + exp(2fc7ri02 n ' > ) + ••• + ^vi^kmifi^-i) p 



n — 1 

(Again, the convergence is almost sure, unless /i is uniform on C.) 



Ck- 



Finally, using arguments analogous to those in the proof of Weyl's equidistribution 



criterion, we shall arrive at our final result. 
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Assume, initially, that \x is not the uniform law on C. For the first task as noted 
above, observe that, by Lemma [2.21 given any e > 0, 

n-l 



3=1 

Now, for any fixed positive integer k, (1 — e ) fc l{ f .(«) e [ 1 _ e qy — (i r j) k < 1, an d so 

j n— 1 j n— 1 



frj n) ) fe < 1. 



3=1 3=1 

Clearly then, a simple squeeze theorem argument gives us 

n-l 



< 5 > ^tB 

3=1 



Now, from Proposition I3.2[ for any positive integer k, 

{y (n) )k + {y (n ))k + + (y W i)fc 



a.s. 

, ► Cfc, 

n — 1 

(r[ n) ) fc exp(2^^ n) ) + (4 n) ) fc exp(2fc7ri(/4 n) ) + ... + (rj n J 1 ) fc exp(2A;7r^g 1 ) a 

n — 1 

Then, applying Lemma 1331 with {x nm = 1 — r^}, we get 

exp(2/c7ri0^ n) ) + exp(2fc7ii02 n) ) + ... + exp(2/c7ri0„_ 1 ) tt . s . 

(6) > c k . 

n — 1 

Now, for the final stage of our proof, 
c fc = E(Z k ), where, Z ~ /i. 
=>• ^ = E(exp(2A;7ri0)) = E(cos(2/c7r6)) + iE(sin(2A;7re)), where, 9 ~ v. 
So, dBJ) gives, 

cos(2fc7r0i n) ) + cos(2fc7r^ n) ) + ... + cos^tt^) a.s. n „ 

>■ E(cos(2A;7rB)), 

n 

sirugwg) + sin(2fc7r^ n) ) + ... + sin^Tn/S) . 

n 

Then, for any trigonometric polynomial, 

r 

q{x) = ao + yj(«j cos(27rjx) + 6j sin(27rjx)), a,, 6j e R, j = 1, 2, r, 

3=1 



Cfc- 
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we have, 



(7) 



n 



E(g(6)). 



Let / be a continuous real- valued function on [0, 1] and fix e > 0. Stone- Weierstrass 
theorem ( [St48] ) states that for any compact Hausdorff space X, and any subalgebra 
S of the space C(X, R) (continuous functions from X to R) that separates points and 
contains at least one non-zero function, S will be dense in C(X, R) with respect to 
uniform norm. In particular, this means that trigonometric polynomials are dense 
in the space of continuous functions on [0, 1] with the uniform norm. Hence, there 
exists a trigonometric polynomial q such that \f — q\ < e. So 



n 



E(/(©)) 



< 



+ 



< 



+ 



"- 1 ,,(4 n) ) 



n 



n 



n 



E(g(0)) 



n 

E-=^(0f h 



n 



-E(?(6)) 



+ |E( ff (6))-E(/(e))| 



+ E| g (e)-/(e)|. 



The first and third terms on the right hand side are each < e while the second term 
goes to almost surely, by ((7j). Hence for any / continuous on [0, 1], 



(8) 



E 



n-l 
J 



n 



E(/(6)), 



and this holds for complex-valued continuous functions as well (which is easily seen 
by comparing the real and imaginary parts). Thus, the joint empirical distribution 
of 0j n \ j = 1, 2, n — 1, converges weakly to is, which means that the joint empirical 



distribution of exp(27T2<; 



1,2, ...,n — 1, converges weakly to /i. This, along 



with Lemma [2.21 gives us the desired result for jj, not uniform on C. 
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Now suppose [x is the uniform law on the unit circle. Then, by @, for any 5 > 0, 



n 



j=l / \ j=l 

S-l + (l-e 



n-l , , ,., x fc 



n- 1 — i' j e[i- 6 ,i]} - (1 -e) fc 
j=i 

We may choose e in such a way that 1 — (1 — e) k < 5/2, and so, the right hand 
side of the above inequality goes to 0. Thus, for any positive integer k, 

n-l 



n - 1 ^ ' 3 



1 ^ 

Note that the above is a slightly weaker version of (151) , since the convergence is 
now in probability, and not almost sure. 

For the rest of the proof, observe that we can follow the same arguments as in the 
non-uniform case, except that the almost sure convergence in each of the statements 
will be replaced by convergence in probability. Thus, we shall arrive at 

5^Ae(/(6)), 

n 

for any continuous function / : [0, 1] — > C. Then, as before, the joint empirical 
distribution of <f>j\j = l,2,...,n — 1, converges weakly to v (which is the uniform 

law on [0, 1]), and so, the joint empirical distribution of exp(27ri(/>^), j = 1, 2, n— 1, 
converges weakly to uniform on C. Lemma [2.21 then gives us the desired result. □ 
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